A new analysis of 14 beta decay: branching ratios and CVC consistency 
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The ground-state Gamow- Teller transition in the decay of 14 is strongly hindered and the electron 
spectrum shape deviates markedly from the allowed shape. A reanalysis of the only available data on 
this spectrum changes the branching ratio assigned to this transition by seven standard deviations: 
our new result is (0.54±0.02)%. The Kurie plot data from two earlier publications are also examined 
and a revision to their published branching ratios is recommended. The required nuclear matrix 
elements are calculated with the shell model and, for the first time, consistency is obtained between 
the Ml matrix element deduced from the analog gamma transition in 14 N and that deduced from 
the slope in the shape-correction function in the beta transition, a requirement of the conserved 
vector current hypothesis. This consistency is only obtained, however, if renormalized rather than 
free-nucleon operators are used in the shell-model calculations. In the mirror decay of 14 C a similar 
situation occurs. Consistency between the 14 C lifetime, the slope of the shape-correction function 
and the Ml matrix element from gamma decay can only be achieved with renormalized operators 
in the shell-model calculation. 

PACS numbers: 23.40.Hc, 21.60.Cs, 21.10.Tg, 27.20.+n 



I. INTRODUCTION 

The nucleus 14 decays predominantly by a Fermi 
transition to the 2.313 MeV 0+ state in 14 N. Weak 
Gamow- Teller branches are evident to the 1 + state at 
3.95 MeV in 14 N (branching ratio, R GT = 0.055%) and 
to the 1 + ground state with Rqt ~ 0.6%. It can be con- 
cluded then that the Fermi transition has a branching 
ratio of ~ 99.3%. Since systematic studies of Fermi su- 
perallowed transitions require their branching ratios be 
known to an accuracy of ±0.1%, the branching ratio of 
the ground-state Gamow- Teller transition for 14 O must 
be determined to within 10% of its central value. 

This Gamow- Teller transition is strongly inhibited. Its 
/f-value is roughly 10 times larger than is typical for 
favoured + — > 1 + transitions. (Even more inhibition 
is evident in the analog 14 C — ► 14 N transition.) The 
inhibition is attributed to accidental cancellation in the 
allowed Gamow- Teller matrix element for this transition 
0. Because the allowed matrix elements are so small, 
the induced terms (particularly "weak magnetism"), as 
well as the relativistic and the second-forbidden terms are 
expected to contribute appreciably to the decay proba- 
bility. As a consequence, many of the usual assumptions 
in the allowed approximation may not be valid. For ex- 
ample, the spectrum shape may deviate markedly from 
the allowed (or statistical) spectrum shape. 

To date there has only been one measurement, by 
Sidhu and Gerhart 0, 0], of the detailed shape of the 
beta spectrum from 14 decay, from which they deter- 
mined the branching ratio of the ground-state Gamow- 
Teller transition with the required precision. It was per- 
formed with an iron-free, beta-ray spectrometer, and was 



published 40 years ago! More recently, calculations by 
Garcia and Brown [j] could not satisfactorily fit the ob- 
served beta spectrum, which led the authors to suggest 
there might be some systematic problem with Sidhu and 
Gerhart's measurement. This conclusion would have a 
serious impact on the branching ratio for the Fermi tran- 
sition. For this reason, we have reanalyzed the data of 
Sidhu and Gerhart, which we obtained from a copy of 
Sidhu's Ph.D. thesis Our conclusion is that the 14 
spectrum shape can be understood, but only if renormal- 
ized operators are used in the shell-model calculations 
of the nuclear matrix elements. Our re-analysis yields a 
ground-state branching ratio of (0.54±0.02)%, compared 
with the originally published result of (0.61 ± 0.01)% 
- a large shift in terms of the uncertainties quoted. 



II. RE-ANALYSIS OF SIDHU-GERHART 
EXPERIMENT 

The probability per unit time for emission of a positron 
whose momentum lies between p and p + dp is 



dN 
dp 



G 2 V 2 

*—* TP * It. 



p\E Q -Ef 



F(Z,E)S(Z,E)Q(p)R(p) 



(1) 
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where G F is the weak interaction coupling constant, V u d 
the up-down matrix element of the CKM matrix, E 
the positron kinetic energy in MeV units, E its max- 
imum value, p the positron momentum, F(Z, E) the 
Fermi function, Z the atomic number of the daughter 
nucleus, S(Z, E) the shape-correction function, Q(p) an 
atomic screening correction, and R(p) a radiative correc- 
tion. (The average value of R(p), when integrated over 
an allowed electron spectrum, is denoted S' R in our pub- 
lications Q on superallowed Fermi transitions.) In this 
work, the shape-correction function includes the nuclear 
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matrix elements 1 . We will write 
27r 3 (n C ) 6 ?iln2 _ — 



G?V2,(m e c»)6 



= K = 6146 s 



(2) 



and let all the instrumental factors be encompassed by 
some function K(j>). For an iron- free beta spectrometer, 
Sidhu assumes the instrumental function is proportional 
to the positron momentum, that is K(p) = up. If dN/dp 
is now interpreted as the number of counts in the beta-ray 
spectrometer per unit time that correspond to positrons 
whose momentum is between p and p + dp, we can write 



dN 
dp 



Ap 3 (E Q - E) 2 F(Z,E)S(Z 1 E)Q(p)R(p), (3) 



where A = K\n.2/K(m e c) 5 and is a constant of unknown 
magnitude. To calibrate the spectrometer and determine 
A, Sidhu made measurements of the lower-energy Fermi 
transition, for which all the factors on the right-hand-side 
of Eq. © are calculable, including the shape-correction 
function S(Z,E). In Fig. 22 of his thesis, Sidhu showed 
the measured Kurie plot for the Fermi transition defined 



1(E) 



I dN/dp 
p 3 F(Z, E) 



A 1 / 2 (E Q -E) 
[S(Z,E)Q(p)R(p)}^ 2 . (4) 



In selecting an energy at which A is to be determined, 
Sidhu J3J writes "we want to avoid being too close to the 
end-point, because of the larger proportion of the high- 
energy positrons (from the Gamow- Teller decay) being 
mixed in this branch near the end point. On the other 
hand, we want to avoid positrons due to contaminants 
such as 15 O. The point at 1.691 MeV should be a good 
compromise." At E — 1.691 MeV, after reducing his re- 
sult by 1.7% to eliminate the estimated contribution of 
15 contaminants, Sidhu quotes 2 



(Eq - Ef 



1.232 ± 0.013 = AS(Z,E)Q(p)R(p). (5) 



Without the benefit of today's computing facilties, Sidhu 
assumed the product S(Z, E)Q(p)R(p) to be equal to 
2, the square of Fermi matrix element. We have 
now computed these factors more exactly and obtained: 
S(Z,E) = 1.00226 x 2.0, Q(p) = 1.00058 and R(p) = 



1 Further, these nuclear matrix elements include the coupling con- 
stants in their definition. For example, the Gamow- Teller matrix 
element, Mqx, includes the axial- vector coupling constant g A . 

2 Actually he quotes a value of 1.237±0.013, but a trivial numerical 
error on page 138 of the thesis indicates that one entry there 
is incorrectly recorded. We have worked backwards from the 
published result Q to deduce what this value needs to be. Note 
the discrepancy is within the stated error. 



TABLE I: Corrected shape-correction functions, C(Z, W), ob- 
tained from Sidhu and Gerhart's data, J(E), via Eq. JHJ. 



E J(E) 
MeV Units: 10~ 



F L /F Q( P ) 



R(p) C(Z, W) 
Units:l(T 4 



2.235 
2.395 
2.544 
2.702 
2.860 
3.018 
3.172 
3.323 
3.482 
3.640 
3.798 



1.738(10) 
1.720(10) 
1.669(10) 
1.652(10) 
1.616(10) 
1.602(10) 
1.553(10) 
1.537(10) 
1.515(10) 
1.500(13) 
1.480(17) 



1.00296 
1.00309 
1.00321 
1.00334 
1.00347 
1.00359 
1.00372 
1.00384 
1.00397 
1.00410 
1.00423 



1.00046 
1.00043 
1.00041 
1.00039 
1.00037 
1.00035 
1.00034 
1.00033 
1.00031 
1.00030 
1.00029 



1.00907 
1.00736 
1.00571 
1.00388 
1.00195 
0.99987 
0.99767 
0.99527 
0.99239 
0.98893 
0.98443 



2.781(16) 
2.755(16) 
2.679(16) 
2.656(16) 
2.603(16) 
2.585(16) 
2.512(16) 
2.492(16) 
2.463(16) 
2.447(21) 
2.424(28) 



0.99523. This product is 0.2% less than Sidhu's choice 
and yields 



A = 0.6173 ±0.0064. 



(6) 



With A now determined, we consider the higher- 
energy lower-intensity Gamow- Teller transition. For this, 
we need to perform a nuclear-structure calculation of 
S(Z, E) and compare with the experimental data. To 
this end it is convenient to introduce the shape-correction 
function defined by Behrens and Biihring 



C(Z,W) = ^^S(Z,E), 



(>) 



where W is the total positron energy in electron rest- 
mass units, W = 1 + E/(m e c 2 ). The factor F/FqLo 
corrects for the fact that Behrens and Biihring compute 
the electron density at the origin and not at the nu- 
clear radius R, as was historically the case. In Fig. 20 
of his thesis, Sidhu plots his experimental data for the 
shape-correction function (corrected for backscattering) 
for eleven positron kinetic energies between 2.2 and 3.8 
MeV. Two lower energy points are disregarded as being 
unreliable. We convert these data to the Behrens and 
Buhring-defined shape-correction function, and further 
apply a screening and radiative correction as follows: 



C(Z, W) 



J(E) F(Z,E) 
AQ(p)R(p) F Q L ■ 



(8) 



Here we have labelled Sidhu's original data as J(E). In 
Table [I] we give J(E) as well as the corrections F/FoLq, 
Q(p), R(p) and the determined values of C(Z, W), which 
are also plotted in Fig. ^ 



III. CALCULATION OF C(Z, W) 

We will use the formalism of Behrens and Biihring [f| 
for the shape-correction function, where it is written in 
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terms of amplitudes MK{k e , k v ) and txik {k e , k v ): 
C(Z,W) = ^k e {M 2 K {kM+m 2 K {kM 

k c k„K 

- \ k '^ M K (k e ,k v )m K (k e , j . (9) 

Here and fik e are beta-decay Coulomb functions, 
which depend on the amplitudes of the electron wave 
functions at the origin, and are defined such that their 
values are of order unity, with corrections of order (aZ) 2 . 
These quantities have been tabulated by Behrens and 
Janecke 0. The factor j kc is defined as [k 2 — (aZ) 2 ] 1 / 2 . 
Here k e and fc„ are the partial-wave expansion labels for 
the electron and neutrino wave functions. In our evalu- 
ations we will keep the lowest two partial waves in each. 
Finally, K is the multipolarity of the transition operators. 
For the ground-state Gamow- Teller decay, + — > 1 + , the 
multipolarity is restricted to K = 1 only. 

Behrens and Biihring [fj give approximate expressions 
for C(Z, W) by expanding the electron wave functions in 
a power series in WR and (cnZ), and the neutrino wave 
function in a power series in p„R, where p v is the neu- 
trino momentum, and R is a typical nuclear size param- 
eter. We have not followed this procedure, but rather we 
have computed Eq. exactly by solving the Dirac equa- 
tion for electron wave functions as described in Ref. @ . 
Our only simplifying assumption is in the evaluation of 
Mj((k e , fc„) and mjf(fc e , k u ), where the nuclear form fac- 
tors of Behrens and Biihring are replaced in the 'impulse 
approximation' by nuclear reduced matrix elements. We 
compute these reduced matrix elements within the nu- 
clear shell model. 

For A = \A nuclei^we use the Op-shell wave functions of 
Cohen and Kurath |9j denoted (8-16)POT. For compar- 
ison purposes we also consider the more recent Op-shell 
part of the Warburton-Brown Hamiltonian El (the in- 
teraction labelled PWBT in Table X of Ref . fig ) . This 
interaction was determined from a least squares fit to 51 
Op-shell binding energies for which the rms deviation of 
the fit was 378 keV. 

These two sets of interactions only incorporate the Op- 
shell. We wanted also to examine the possible effects 
of s<i-shell contributions. Close to major shell closures 
the choice of a model space and effective interaction can 
be problematic if one wants to go beyond simple single- 
major-shell configurations. For example, in the A = 14 
spectrum the lowest-energy states are predominantly two 
holes outside a closed 16 core, \2h), but lying low in 
the spectrum are 'intruder' states with configurations in- 
volving four holes and two particles, |4ft-2p). Mixing be- 
tween these configurations must occur, and to obtain the 
degree of mixing with the shell model is difficult. Shell- 
model calculations that attempt to mix |2ft) and |4ft-2p) 
configurations encounter what has been called [l(j the 
"nftw catastrophe". The presence of |4ft-2p) configura- 
tions depresses the \2h) states, opening up a large energy 
gap between the \2h) and |4ft-2p) states. This would 



be corrected somewhat if the model calculation included 
|6ft-4p) states as well, since the role of the |6ft-4p) states 
is to depress the |4ft-2p) states. Thus when truncating 
the model space to include only \2h) and |4ft-2p) states, 
the depression driven by the |6ft-4p) states on the |4ft-2p) 
states is absent. In an attempt to circumvent this catas- 
trophe we will use quite weak cross-shell interactions and 
examine the sensitivity of our results to the strength of 
the cross-shell interaction. 

We adopt the following method for incorporating sd- 
shell effects in the mass- 14 system: We use the Cohen- 
Kurath interaction [9j for p-shell interactions, the USD 
llj for s,d-shell interactions and the Millener-Kurath 
|13 | interaction for the cross-shell matrix elements. The 
Millener-Kurath interaction was designed to reproduce 
the unnatural-parity states in p-shell nuclei, such as the 
negative-parity states in A = 14 that involve just one 
particle in the s, d-shell. It wasn't designed to give the 
mixing between \2h) and |4ft-2p) configurations. Nev- 
ertheless we will use the (2h\V\4h-2p) matrix elements 
given by the Millener-Kurath interaction and multiply 
the matrix elements by a factor, /, that ranges from 0.0 
to 0.6. When / = 0.0, there is no mixing between the 
1 2ft) and |4ft-2p) configurations, and when / = 0.6 the 
ground-state wave function is approximately 74% |2ft) 
and 26% |4ft-2p). We will quote results using /=0.3, and 
denote this interaction as MK. We have examined the 
sensitivity of our results to a variation of / and found 
that the spread of the different results is within the as- 
signed errors. 

We are interested in further refining the wave function 
for the 1 + T = state in A = 14. As was noted by Garcia 
and Brown p|, who suggested this procedure, the GT 
transition strength to the lowest 1 + state is very small 
compared to the strength of the transition to the second 
1 + state at 3.95 MeV excitation energy. Thus, any small 
mixing between these two 1 + states in the model will 
have a large effect on the weak transition rate. To make 
use of this fact, we can write the wave function for the 
lowest 1 + state in 14 N as 

|l+low)=a|l+(l))+/3|l+(2)) (10) 

with a 2 + (3 2 — 1. Here (1) and (2) refer to the first and 
second model states obtained with either the CKPOT, 
PWBT or MK effective Hamiltonians. In fitting the beta- 
decay data, it turns out that we need a negative sign for 
(3 with the CKPOT interaction and a positive sign with 
the PWBT and MK interactions. 

Our strategy is then to adjust a to minimize the \ 2 
between the calculated C (Z, W) and the corrected exper- 
imental shape-correction function given in Table [I] The 
Gamow- Teller matrix element, Mqt, is particularly sen- 
sitive to small variations in a, and consequently is rather 
precisely determined mainly from the fit to the absolute 
magnitude to the shape-correction function. We have 
fitted our calculated C(Z, W) to the expression 

C(Z, W) = |M GT | 2 fc(l + aW + mjib/W + cW 2 ) (11) 
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TABLE II: The value of the wave function amplitude, a, Eq. IIUH that gave the best fit to the experimental shape-correction 
function for 14 0, and the best fit to the decay half-life for 14 C. Also given are the Gamow- Teller and Ml matrix elements for 
the + — > 1 + transition, the parameters k and a of the shape-correction function, C(Z,W), and the ground-state branching 
ratio, Rgt- The \ 2 P er degree of freedom, \ 2 /v for the fit to the 14 shape-correction function is recorded. Free-nucleon 
operators were used for the shell-model calculations. 



Model 



X 2 /u M G i 



a 

MeV" 



Rgt(%) 



14 decay: 

CK 0.98447 
PWBT 0.99805 
MK 0.98482 



2.0 0.01400 
3.5 0.01392 
1.5 0.01406 



-0.533 
-0.569 
-0.515 



1.999 
2.068 
1.958 



-0.133 
-0.139 
-0.130 



0.566 
0.570 
0.563 



14 C decay: 

CK 0.98560 
PWBT 0.99760 
MK 0.98595 



-0.00465 
-0.00469 
-0.00458 



0.503 
0.539 
0.484 



0.345 
0.346 
0.354 



-0.647 
-0.683 
-0.624 



Expt: 
Expt: 



MffJ from T 7 in 14 N 
4 C slope a 



0.312(7) 



-0.45(4) 



as the approximate expressions of Behrens and Buhring 
can be cast in this form. In Table ILTI we give the first two 
parameters, k and a, as determined by least-squares fit- 
ting for each set of shell-model interactions. The param- 
eter a, which is called the slope of the shape-correction 
function, can be expressed in approximate form 



M 





Ml 



3M M GT 



as" 



(12) 



where M is the nucleon mass and Mffj is the Ml ma- 
trix element. Clearly this slope is dominated by an in- 
terference between the Gamow- Teller and Ml matrix ele- 
ments. This is of considerable - indeed historical - impor- 
tance 13]. A measurement of M^ I1 obtained from beta 
decay can be compared with the corresponding matrix 
element, Mj 41 , obtained from the electromagnetic tran- 
sition between the isobaric analog state and the ground 
state of 14 N. The conserved vector current hypothesis as 
enunciated by Gell-Mann [l3j proposes that the vector 
current of the weak interaction is just a rotation in isospin 
of the vector current of the electromagnetic interaction; 
in this case we would expect M^ n = — M^/%/2. (The 
V2 factor originates in the isospin Clebsch-Gordan coef- 
ficient for a charge-changing reaction, such as beta-decay, 
differing by \[2 from the isospin Clebsch-Gordan coeffi- 
cient for a charge-conserving reaction, such as gamma- 



3 This approximation is quite poor. We only display it to show 
the important effect of the Ml matrix element. 



decay.) The analogous electromagnetic transition in 14 N 
is the 0+ T = 1 -> 1+ T = decay of the 2.313 
MeV state, which has a measured gamma width of 
T 7 = (6.7 ± 0.3) x 1CT 9 MeV. The corresponding Ml 
matrix element is 



I M M1 1 expt 



3I\ 



= 0.442 ±0.010. (13) 



Here E 1 is the photon energy in MeV, E 1 = 2.313 MeV, 
and /j,n the nuclear magneton unit, fi^ — 0.1262 MeV 1 / 2 
fm 3 / 2 . The matrix element is dimensionless. To conform 
with CVC, the beta-decay Ml matrix element should be 
given by 



|M 



I 

Ml I 



\Ml [1 \/V2 = 0.312 ±0.007, 



(14) 



if we neglect charge-symmetry breaking corrections. 
Garcia and Brown jj] have studied this issue and con- 
cluded that "one should not expect very large charge- 
symmetry-breaking effects in the A = 14 system" . 

It is evident in Table ILTI that for values of a that give a 
good fit to the experimental beta-decay spectrum shape 
the value of the Ml matrix element differs from Eq. i|14|) 
by about a factor of 1.7. This disappointing result is also 
consistent with the conclusions of Garcia and Brown Q : 
with Op-shell wave functions it is not possible to fit the 
5(M1;0 + — > 1 + ) radiative decay simultaneously with 
the beta-decay measurements of Sidhu and Gerhart [3-E1 

We have also examined the mirror decay of 14 C, whose 
lifetime is known ^^] to be ti/ 2 = 5700(30) yr. Our strat- 
egy for this decay is to adjust the wavefunction amplitude 
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TABLE III: Same as Table ITT1 except that renormalized operators are used rather than free-nucleon operators in the shell-model 
calculations. 



Model 



X 2 /u Mai 



a Rgt{%) 
MeV -1 



14 O decay: 

CK 0.98428 5.1 0.01490 -0.313 1.552 -0.099 0.540 

PWBT 0.99812 3.3 0.01480 -0.346 1.605 -0.105 0.544 

MK 0.98463 6.0 0.01494 -0.300 1.526 -0.095 0.538 



14 C decay: 

CK 0.98556 
PWBT 0.99762 
MK 0.98590 



-0.00338 
-0.00343 
-0.00333 



0.281 
0.314 
0.267 



0.613 
0.607 
0.629 



-0.421 
-0.459 
-0.406 



Expt: |Mf n | from T 7 in 14 N 
Expt: 14 C slope a 



0.312(7) 



-0.45(4) 



a to fit this lifetime exactly. Then the shell-model cal- 
culations give a prediction for the shape-correction func- 
tion, C(Z, W), whose slope parameter a can be compared 
with the experimentally determined value of Wietfeldt et 
al. [l^. The results are listed in the bottom half of Ta- 
ble [D] Again we have a disappointing result: the calcu- 
lated slope parameter of a ~ —0.65 MeV -1 differs signifi- 
cantly from the experimental value 4 of —0.45(4) MeV" 1 . 

A. Renormalized operators 

We carried out the shell-model calculations just dis- 
cussed using operators derived in the impulse approxi- 
mation with coupling constants appropriate for free nu- 
cleons. In finite nuclei one expects corrections to this 
scheme coming from two sources: firstly, the shell-model 
calculation is carried out in a truncated model space, 
which can be corrected in a perturbation expansion, and 
secondly the nucleons in the nucleus are interacting via 
the exchange of mesons and these mesons can influence 
the electromagnetic and weak interactions in nuclei. The 
two corrections are called core polarization and meson- 
exchange currents respectively. These phenomena are re- 
sponsible for the quenching of the Gamow- Teller matrix 
element in finite nuclei. We define 

ffA.cff = ,9a + 5g A , (15) 



4 We take the slope parameter from Table II of llfl using the fitted 
value for the electron spectrum in the 100 to 160 keV energy 
range. In a wider energy range, 50 — 160 keV, the authors quote 
a slope parameter of —0.32 MeV -1 . 



where g A is the free-nucleon value of the axial- vector cou- 
pling constant, g A ~ 1.27, and Sg A the correction to it. 
We fix the value of Sg A by considering the beta decay 
of 15 to its ground-state mirror in 15 N, whose experi- 
mentally determined Gamow- Teller matrix element |l7j 
shows a reduction of 13.2(7)% over that calculated with 
the free-nucleon coupling constant for a configuration 
that consists of a single Op-hole in a closed 16 O core. 
To fit the experimental value requires 

6g A = -0.165, (16) 

and we will adopt this value. 

We also need to understand how the isovector Ml op- 
erator is renormalized by core polarization and meson- 
exchange currents. Both these corrections have been 
evaluated in an ab initio calculation for closed-shell- 
plus-or-minus-one nucleon configuration by Towner and 
Khanna [lg. Their results for the A — 15 case of a Op 
hole in an 16 O core are expressed in terms of an effective 
Ml operator defined as 

(Ml) cff = gi%L + gi%S + gi%[Y 2 , S] , (17) 

where gx^a — fi4 + $9x\ with g x being the free-nucleon 
coupling constant, Sg x the calculated correction to it and 
X = L, S or P. The free-nucleon values are g^ =0.5, 
gi 1 ^ = 4.706 and g^ = 0. Towner and Khanna cal- 
culated the isovector combinations of the magnetic mo- 
ments of the 15 and 15 N ground states and obtained 
an 8.9% enhancement in the free-nucleon isovector mag- 
netic moment compared to an experimental enhancement 
of 11.1% - a clear success for an ab initio calculation. We, 
therefore, adopt their calculated values of 

SgP = 0.076 SgP = -0.22 Sg^ = 0.96. (18) 



6 



C(Z,W) 

XlO 4 
2.8 



2.7 



2.4 



1 1 1 1 1 1 

\. free-nucleon operators 


1 1 


" X\ 




- \ 




renormalized operators 
i i i i i i 


X 



2.3 2.5 2.7 2.9 3.1 3.3 3.5 3.7 
Positron kinetic energy (MeV) 



TABLE IV: Statistical rate function /, Eq. for the 

Gamow- Teller transitions in 14 and 14 C decays. 







/ 


Model 


14q 


14 C 



free-nucleon operators 

2473.1 1.5552 x 10" 



CK 
PWBT 

MK 



2522.9 1.5284 x 10 
2441.8 1.6040 x 10 



renormalized operators 

CK 2086.4 2.9314 x 10 
PWBT 2128.1 2.8437 x 10 
MK 2067.5 3.0222 x 10 



FIG. 1: The shape-correction function calculated with PWBT 
wave functions compared with the experimental data of Sidhu 
and Gerhart Q. 



1633.6 6.0766 x 10" 



With the effective operators thus selected, we repeated 
the strategy of adjusting the wave function amplitude 
a to minimize the \ 2 between the calculated C(Z, W) 
and the corrected experimental shape-correction function 
given m Table H The results of the fit are given in Ta- 
r/3 



blc lIIII where it is observed that the M^ n matrix element 
is considerably reduced from that in Tableland is now 
comparable to that deduced from gamma decay. This is 
a major success and a strong endorsement for the use of 
renormalized operators. However, it is not all good news. 
For two of the wave function choices, CK and MK, the 
quality of the fit to the experimental shape-correction 
function is inferior, the \ 2 P er degree of freedom being 
two to four times larger. Only for the PWBT interaction 
is the quality of the fit comparable. We show the PWBT 
fit in Fig. ^ It is clear that the renormalized opera- 
tors lead to a smaller slope in the shape-correction func- 
tion. This is the expected result since the slope is gov- 
erned by the matrix-element ratio M^/Mqt, as given 
in Eq. (Q2J. 

For the mirror decay in C, our strategy was to fix 
the wavefunction amplitude a to reproduce the known 
lifetime, and use these wave functions to compute the 
shape-correction function C(Z,W), giving a prediction 
for the slope parameter, a. From Table ITTT1 it is clear that 
the use of renormalized operators is very successful: the 
calculated slope parameter of a ~ —0.43 MeV -1 agrees 
perfectly with the experimental result of Rcf. 16] . 



IV. 



THE 14 GAMOW-TELLER BRANCHING 
RATIO 



Now that a shape-correction function C(Z,W) has 
been obtained for 14 decay that agrees reasonably well 



with the data of Sidhu and Gerhart 0, Q (see Fig. ^| we 
integrate this function over the entire beta spectrum and 
compare with the analogous Fermi transition to obtain 
the Gamow- Teller branching ratio, Rgt- 



R 

Rgt 



t_F _ 
tGT 

R 
l + R' 



N GT /gt(1 + ^)gt|M gt | 2 



N F / F (1 + ^) F |M F | 2 (1 



s c y 

(19) 



where tp and tar are the partial half-lives of the Fermi 
and Gamow- Teller transitions respectively, and Np and 
Nqt are their integrated count rates, N = J (dN/dp)dp. 
For the Gamow- Teller transition, the statistical rate func- 
tion is defined as 



1 



\M G1 



pW(W - W) 2 F L C{Z, W)Q(p)dW. 

(20) 

Note that, since C[Z, W) includes the nuclear matrix 
elements, we have divided by lAfc^l 2 to conform with 
the normal definition of /. The calculated values for the 
Gamow- Teller transition are listed in Table IIVI for the 
various wave function selections. 

The other factors in Eq. (|19f) were evaluated as fol- 
lows: For the Fermi transition we calculate fp — 42.772 
(see Ref. Q) and |M F P = 2.0. Next, using the meth- 
ods described in Ref. [a] we obtain (S' r )gt = 1.294%, 
(S' R ) F = 1.520% and 5 C = 0.57%. With these values we 
obtain the branching ratios listed in the last column of 
Tables ITT1 and IIIII We adopt the result with renormalized 
operators as our central value and assign an error that 
is half the spread between the renormalized and free- 
nucleon operators, obtaining Rgt = (0.540 ± 0.015)% . 
However, there still remains a 1% normalization uncer- 
tainty in the value of the calibration constant A, Eq. JfJJl. 
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We add this uncertainty linearly to obtain as our final 
branching ratio a value of 

Rgt — (0.54 it 0.02)%. (21) 

This result differs from Sidhu and Gerhart's published 
result of Rqt = (0.61±0.01)% and has a larger uncer- 
tainty. The Gamow- Teller matrix clement from the fit is 
M GT = 0.0149±0.0005 compared to Sidhu and Gerhart's 
published result of M GT = 0.0164 ± 0.0004. 

V. IMPACT ON SUPERALLOWED BRANCH 
A. Two earlier measurements 

There are two earlier measurements of the 14 O branch- 
ing ratio: Sherr et al. 0] obtained Rqt = (0.6 ±0.1)%, 
while Frick et al. HJ obtained (0.65 ± 0.05)%. In both 
cases Kurie plots were constructed from raw data for both 
the Fermi and Gamow- Teller transitions and the required 
branching ratio was obtained from their ratio: 



R 



If X% 



Rgt 



R 



1 + i?' 



(22) 



where Xf, Xqt are the ratio of the Kurie plot data to 
the allowed approximation (Wq — W), viz. 



1=1 



K(Wi) 



K(W) 



I dN/dp 
p 2 F(Z, W) ' 



(23) 



Here Wi are the values of W for which experimental data 
have been obtained, and rij is the number of such data. 
The unknown normalization of the Kurie plots cancels in 
the ratio. Further, / is the integrated electron spectrum, 
which in the allowed approximation is 



/stat — 



pW(W - W) 2 F(Z, W)Q{p)dW. (24) 



Sherr et al. [l!j only published their Gamow- Teller Kurie 
plot, so it is not possible to reanalyze their result. Frick 
et al. [20j, on the other hand, published both their Fermi 
and Gamow- Teller Kurie plots, so we have reanalyzed 
them according to Eq. (|2*2)> obtaining Rqt — (0.64 ± 
0.03)% agreeing satisfactorily with the published result. 
(Our uncertainty only includes the statistical uncertainty 
in the fit of the allowed shape to the Kurie data, and the 
uncertainty in the end point energy, Wo, as known in 
1963). 

The original analyses of both experiments were based 
on the allowed approximation. However, as we have 
discussed at length, the 14 Gamow- Teller transition is 
strongly hindered and the allowed approximation not suf- 
ficient. Thus, we have reanalyzed the Frick et al. |2£j 
Kurie plot data using for X 



X 



1 th 

n. ^ — ' 



K{Wi)M t 



GT 



n * ti ( w o - w i) Vc(^W) 



(25) 



and using for the integrated spectrum, /, the exact ex- 
pression given in Eq. (|20[l . The shape-correction func- 
tion, C(Z,W), is evaluated with p-shell wave functions, 
PWBT, and renormalized operators. The result is an in- 
crease in the branching ratio to Rgt — (0.73 ± 0.03)%. 
We checked our procedure by performing the same anal- 
ysis on Sidhu's Kurie plots, as published in his thesis 
[3j. In the allowed approximation we obtained Rgt = 
(0.45 ± 0.01)%, while with the exact expressions we ob- 
tained Rgt = (0.53±0.01)% in agreement with our more 
accurate analysis of the shape-correction functions. The 
conclusion is clear. For the Gamow- Teller branching ratio 
in 14 O, determinations based on an allowed approxima- 
tion analysis of Kurie plots have to be increased by about 
14%. Unfortunately this places the earlier data in conflict 
with the more accurate Sidhu-Gerhart result |2j,|3jj. This 
is not a question of the method of analysis: the raw data 
are in conflict. If one considers a ratio of ratios, compar- 
ing the ratio of Gamow- Teller to Fermi Kurie plots for 
Sidhu Q and Frick et al. [2(j , discrepacies of order 20% 
are evident. This, alone, leads to a 40% difference in the 
deduced branching ratios. If a modern day experiment 
were to be mounted, this would be one discrepancy that 
could quickly be resolved. 



The Fermi branch and corrected Tt value 



A survey of all the data on superallowed + — > + 
Fermi decay has recently been published by Hardy and 
Towner 8]. For the 14 Fermi branching ratio, the value 
of 99.334(10)% is given there based on the ground-state 
Gamow- Teller branching ratio obtained from Sidhu and 
Gerhart Q averaged with the two older and less precise 
results from 0] and [2(|, and a second Gamow- Teller 
branching ratio to the 3.95 MeV state of 0.0545(19)%. If 
we now replace the Sidhu-Gerhart value with the result 
from Eq. I|21f) . and increase the branching for the two 
older data by 14% as discussed in Sect. IV Al but leave 
their error assignments at their published values, then 
we obtain a very conservative estimate of the ground- 
state Gamow- Teller branching ratio of (0.57 ± 0.06)%. 
The uncertainty here has been scaled by 2.6 acoording 
to our usual prescription because of the incompatility 
of the Sidhu-Gerhart result with the two older measure- 
ments. The Fermi branching ratio is now increased to 
99.376(65)%. The impact of this is to lower the cor- 
rected Tt value from 3071.9(26)s to 3070.7(32)s. This 
value still leaves the 14 datum consistent with the other 
eleven precision measured Tt values, although it is now 
on the low side of the average. The slight shift in 14 
is well within the stated errors in the survey and has 
a negligible impact on the physics conclusions obtained 
there. 
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VI. A NOTE ON ft VALUES FOR 
GAMOW-TELLER TRANSITIONS 

It is traditional to characterize an allowed Gamow- 
Teller transition by its log ft value, where the statisti- 
cal rate function / used in this application is devoid of 
any nuclear-structure factors and is defined in Eq. (|24|l . 
When the Gamow- Teller matrix element is large, of or- 
der unity, then the shape-correction function C(Z, W) 
is nearly independent of energy and C(Z, W)/\Mgt\ 2 is 
close to unity. Under these conditions there is little dif- 
ference between the exactly defined / of Eq. (|20[1 and the 
traditional expression in Eq. (|24p. But for the Gamow- 
Teller transitions in 14 and 14 C decays, |Mgt| is very 
small and the shape-correction function has a significant 
effect. For these transitions there is a large difference 
between the exact / and / s t a t as shown in Table llVl 

For the exactly-defined /, the transition ft- value 
equals a constant divided by the square of the Gamow- 
Teller matrix element 



27r 3 ln2(?ic) 6 ^ 1 6146 

^K^) 5 |M GT |2 = ]iW(3Tp 



(26) 



Note that in our notation |Mgt| includes the axial- vector 
coupling constant, g A . For precision work, the lifetime t 
should be adjusted for radiative corrections. Eq. I|26l) is 
frequently used to deduce the Gamow- Teller matrix ele- 
ment from a published ft value for which / s t a t has been 
used for /. For example, the National Nuclear Data Cen- 
ter gives the log ft value for 14 C decay as 9.040(3) 
and for 14 decay as 7.279(8). To deduce \M GT \ from 
Eq. (|26|l with these values would be incorrect. The log ft 
values for the exactly defined / are 8.72(2) and 7.44(1) 
respectively, with larger error bars because of the uncer- 
tainty from nuclear structure. 



VII. CONCLUSIONS 

We began this work disturbed by the statement from 
Garcia and Brown 4] that it was not possible to fit 



the B(M1;0 + — > 1+) radiative decay in 14 N simulta- 
neously with the 14 beta decay measurements of Sidhu 
and Gerhart Q , an apparent violation of the conserved 
vector current hypothesis. A problem with the Sidhu- 
Gerhart work "could significantly change the conclusions 
extracted from + — > + transitions regarding universal- 
ity and unitarity" , they wrote. We initially reanalyzed 
the Sidhu-Gerhart experiment and came to the same con- 
clusion as Garcia and Brown, but we then discovered that 
by using renormalized operators in the shell-model calcu- 
lation we could achieve much greater consistency between 
the requirements of CVC and the measurements of Sidhu 
and Gerhart. This observation was reinforced when we 
examined the mirror 14 C decay. There it was only possi- 
ble to fit the known lifetime and the slope parameter in 
the shape-correction function [ljj while remaining con- 
sistent with the requirements of CVC, if renormalized 
operators were used. This is our principal physics con- 
clusion: the use of renormalized operators is mandatory. 



A second important outcome relates to the superal- 
lowed beta decay of 14 0. With the reanalysis of the 
Sidhu-Gerhart experiment the recommended value from 
this experiment for the Gamow- Teller branching ra- 
tio from 14 to the ground state of 14 N is Rqt = 
(0.54 ± 0.02)% compared to the published value of 
(0.61±0.01)%. This result, when combined with updated 
older measurements, revises the recommended branch- 
ing ratio for the Fermi transition from 99.334(10)% to 
99.376(65)% and shifts the corrected Tt value for the 
Fermi branch downwards by 1.2s. This shift is within 
the stated uncertainty of the Tt value given in the sur- 
vey of Hardy and Towner p| and does not alter any of 
the conclusions reached there. 
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